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Abstract 

Phenomenological Richardson potential has built in asymptotic freedom (AF in short) 
and confinement, with only one parameter A in the potential. But it is known that 
the scales of AF and confinement are not the same. In the present work a relativistic 
mean field calculation for baryons is tried out with two parameters A and A' for AF and 
confinement respectively . 

To test the two parameter potential we calculate the energies and the magnetic mo- 
ments, of the triple u - quark system (A ++ ) and the triple s - quark system (Q~ ) and 
found good values for A = 100 MeV and A' = 350 MeV. So we believe that the 
modified Richardson potential should have AF scale A = 100 MeV and the confinement 
scale A' = 350 MeV. 
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1 Introduction 



't Hooft suggested that the inverse of the number of colors N c could be used as an expansion 
parameter in the otherwise parameter free QCD theory [1]. By the end of 90's, properties 
of large N c baryons have been extensively studied by algebraic methods for spin and isospin 
symmetry. It has become possible to make a unified view on the various effective theories such 
as the Skyrme model, the non-relativistic quark model and the chiral bag model [2]. Witten 
[3] suggested that for large N c baryons a mean field description could be obtained using a 
phenomenological interquark potential tested in the meson sector 1 

Indeed self consistent baryon mass calculation is feasible with success in the mean field level 
[4] using the Richardson potential [5] as an interquark interaction. Richardson potential takes 
care of the two features of the qq force, asymptotic freedom and confinement as given below : 

where — J ^" 1 is the color contribution, Nf is the number of flavors, taken to be three, and 

f (t) - 1 - 4 r - exp (- qt) (2) 

/W - Ji q [ln(<? -1)] 2 + tt 2 [Z) 

A is a parameter whose value was originally chosen as 400 MeV as for small q 2 the potential 
reduces to a linear confinement, and the linear confinement string tension from lattice calcula- 
tion is about that value. However, the asymptotic freedom part also has the same A. In other 
words, both the confinement and asymptotic freedom scales are chosen to be the same in the 
potential equation (1) even though it is known that the asymptotic freedom scale should be 
around 100 MeV. 

Same potential has been included in the relativistic HF calculation of strange quark matter 
to form compact stars. These stars are more compact than the conventional neutron stars and 
fit into the Bodmer- Witten hypothesis for the existence of strange quark matter. They lead to 
systems more compact than those derived assuming the bag model type of quark confinement 
[6, 7, 8] and leads to various observable phenomenon [9, 10, 11, 12, 13, 14, 15, 16, 17, 18]. But 
A parameter had to be within 100 — 200 MeV. This is consistent with asymptotic freedom scale 
but for confinement scale, it is low. 

Since the scales of these two phenomena, need not be the same as in the original potential., 
we incorporate them explicitly as A' and A: 



1 We agree with the referee that a relativistic Faddeev calculation would be more appropriate for a 3 quark 
system but do not attempt that. This is because our goal is to find an effective interaction from finite system 
for use in strange star calculations. 
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A' corresponds to the confinement part and A corresponds to the asymptotic freedom part. 

Recently the magnetic moment of A ++ was derived from experimental data [19] and it was 
shown in [20] that many models of hadrons cannot fit its value satisfactorily. We re-visit the 
relativistic mean field calculation [4] with two parameter Richardson potential (3). Elsewhere 
the star calculation is repeated with the modified potential ([21]). 

In this work an extensive search is performed to find the magnetic moment and mass of 
the A ++ and Q~ . It is found that for the confinement part a A' of about 340-350 MeV and 
for asymptotic part a A of 100 MeV reproduces both the magnetic moment and the energy 
satisfactorily. 

This relatively small value of asymptotic freedom parameter (A) is in agreement with Shif- 
man et al [22] (70 MeV - 100 MeV). 

Interestingly, in a recent paper, Radzhabov and Volkov [23] proposed that A' should be 340 
MeV. 



2 Details of calculation 

Assuming the quarks occupy the same orbital, the mean field w av (r) is given by 



N 2 — 1 r 



(4) 



where — N ^ T 1 is the color contribution in the color singlet state, the Fock term. The color 



2N C 



factor for the Hartree term is zero. The energy E HF is given by 



E HF = N c [e - -< (p jm \uj av \(f) jm 



(5) 



where e is the single particle energy obtained by solving 



[a.p + f3m + u; av (r)]<p jm (r) = e<p jm (r) (6) 
where a's and (5 are the usual Dirac matrices. For quarks in the lowest(lsi/2) orbital, one 



may write 



<^ 1/2 ( r ) 



47T 



iG(r)Xm 
a.rF(r)Xr 



(7) 



where \m is the Pauli spinor and the eqn(5) yields the system of coupled differential equa- 
tions: 



dF 
dr 



dG , 

— - {m-uj av + e)F 



+ (-) F + (e - u av - m)G = 



(8) 
(9) 
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The above equations are solved self-consistently by iteration since u av depends on Dirac 
large and small components G(r) and F(r). 

In eqns (8) and (9), the self consistent single particle confining potential is a vector one, 
leading to instability in the solution. The same problem was encountered by Crater and Van 
Alstine [24] who suggested a prescription of taking a half- vector half-scalar form for the linear 
part. This choice also leads to a cancellation of spin-orbit effects at long range. So we add 
the confining part of the two-body potential equally to the energy and mass (1/2 vector, 1/2 
scalar potential). 

The vector and scaalar potentials are respectively: 



V vec (r 12 ) = 



Nt 



6"7T 



r A' 2 r 12 



2N C 33-2iV7 

T2 



/(Ari 2 ) 



y scalar {f 12) — 



TV 2 -! 



Gn A /2 r 12 



(10) 



(11) 



2N C 33 - 2N f 2 

The vector potential K, ec (r 12 ) is used in the expression of oo av {f) [equation 4] and the scalar 
potential V sca i ar (r 12 ) is added with the mass terms in the coupled differential equations. 
The magnetic moment \i is given by (see [25]):- 



\i = ~e B - I G(r)F(r)r 3 dr 
3 Jo 



(12) 



where is the charge of the baryon. 



It is convenient to express the Dirac components G(r) and F(r) as sum of oscillators : 

G(r) = £C n i?„o (13) 

n 

F(r) = Y, D mRmi (14) 

m 

where C's and D's are coefficients. The expansion reduces the differential equation to an 
eigenvalue problem. The solution is to be found self consistently by diagonalizing the matrix 
and putting back the coefficients till convergence is reached. In general R n i(r) is given by : 



Rni(r) 



j^r l exp(—-r 2 )lJn 2 (r 2 ) 



\T(n + l + l) 



(15) 



i+- 

where L n 2 (r 2 ) are Associate Laguerre polynomials. 

In the calculation, r is replaced by r/b in the expression of Vn 2 (r 2 ) where b is the oscillator 
length which need not to be the same for G(r) and F(r) and they are denoted by b and b' 
respectively. 
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Starting with a trial set of C's and D's , uj av is calculated, a matrix is constructed and 
diagonalized. The new coefficients C n and D n can be read off from the eigenvector with the 
corresponding eigenvalue e. The wave functions are put back in the equations and the process 
is continued till self consistency is reached. 

But the problem is that the RHF solutions violate translational invariance, since they are 
formed by single-particle wave functions derived by an average potential ui av . As a consequence, 
the centre-of-mass momentum is not well-defined in RHF solutions and this entails a spurious 
contribution from the centre-of-mass [CM] kinetic energy to the total energy. Since the relative 
importance of this effect increases as the number of particles decreases, it is important that it 
should be corrected for systems formed of few particles. This can be done by extending to the 
RHF equations the Peierls-Yoccoz procedure of nuclear physics. The spurious contribution is 
denoted by Tqm and the baryon mass [M] has to be compared with the difference Ehf - Tqm- 
Here this spurious contribution has also been calculated and subtracted from E HF to estimate 
the correct energy of the baryon. 

We have also estimated the values of r.m.s. radius r av and checked if the wave functions are 
normalized or not. The normalization factor is denoted by N. 

r av is given by: 



r av = Jj o {G{rf + F{rf)r^dr (16) 

N is given by: 

N = £ maX (G(rf + F(rf)r 2 dr (17) 
where r max is the upper limit of integration and is taken as 5.0 fm to make N ~ 1. 

3 Results for A++ 

We start with the simplest system, the totally symmetric spin and isospin triple u-quark state. 
Here m u is taken as 4 MeV. We checked the convergence in choosing the number of the oscil- 
lators. In figure [1] we see that the change in Ehf is 32 MeV when we increase the matrix 
dimension from 5 X 5 to 7 X 7, but only 16 MeV from 7 X 7 to 9 X 9. It is interesting to note 
that the magnetic moment remains almost same when we increase the dimension from 7X7 
to 9 X 9. In table 1, it is shown that E H f is almost independent of oscillator parameter b and 
b' and the mass after CM correction is 1171 MeV. 

Now, we vary A and A' independently. We give details of A and A' variational results in 
table 2 and table 3 for matrix dimension 7 X 7 as an example. We find that Ehf varies widely 
with A' as expected since A' is the confining parameter. With A' = 350 MeV and A = 100 MeV 
( table [2] ) A mass = 1250 MeV and its magnetic moment 5.77 magneton. Our final result 
is given in table [4] for 9 X 9 matrix. We get the A mass to be 1224 MeV and its magnetic 
moment jj, = 6.15 magneton. 
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Table 1: Variation of Hartree Fock energy, centre of mass kinetic energy, mass and magnetic 
moment of A ++ with oscillator parameter b where A = A' = 350 MeV; for 9 X 9 matrix. 



b' 


b 


E HF (MeV) 


T CM (MeV) 


M(MeV) 


tL 


row (fin) 


N 


0.60 


0.78 


1310 


139 


1171 


6.75 


1.30 


0.988152 


0.60 


0.80 


1307 


138 


1169 


6.69 


1.31 


0.990868 


0.60 


0.82 


1303 


138 


1165 


6.56 


1.32 


0.995152 


0.60 


0.84 


1301 


130 


1171 


6.42 


1.31 


0.998549 


0.60 


0.86 


1299 


130 


1169 


6.30 


1.31 


1.000855 


0.60 


0.88 


1298 


127 


1171 


6.23 


1.30 


1.002296 



1390 I 1 1 1 1 1 

1380 - O o o 

1370 - 

1360 - 
|l350 - 
|l340 - 

U1330 + + + + + 
1320 - 

1310 □ n 

1300 - D □ 

1290 I 1 1 1 1 1 

0.76 0.78 0.8 0.82 0.84 0.86 
Oscillator Parameter (b) 



Figure 1: Variation of Hartree- Fock energy Ehf with Oscillator Parameter b for 5 by 5 , 7 by 7 and 9 by 9 
matrices . Here o stands for 5 by 5 matrices, + stands for 7 by 7 matrices, □ stands for 9 by 9 matrices . 



Table 2: Variation of Hartree Fock energy, centre of mass kinetic energy, mass and magnetic 
moment of A ++ with asymptotic parameter A where A' is 350 MeV; b' is 0.60, b is 0.84 using 
seven by seven matrix. 



A (MeV) 


E HF (MeV) 


7cm (MeV) 


M(MeV) 


LL 


»"av(fm) 


TV 


100 


1391 


141 


1250 


5.77 


1.11 


1.000001 


250 


1340 


141 


1199 


5.95 


1.13 


1.000001 


300 


1333 


140 


1193 


5.99 


1.14 


1.000000 


325 


1330 


140 


1190 


6.00 


1.14 


1.000000 


375 


1325 


140 


1185 


6.02 


1.14 


1.000000 
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Table 3: Variation of Hartree Fock energy, centre of mass kinetic energy, mass and magnetic 
moment of A ++ with confinement parameter A' where A is 350 MeV; b' is 0.60, b is 0.84 using 
seven by seven matrix. 



A'(MeV) 


E HF (MeV) 


T CM (MeV) 


M(MeV) 


it 

MO 


r a „(fm) 


N 


250 


1003 


140 


863 


6.03 


1.39 


0.999999 


300 


1171 


139 


1032 


6.09 


1.22 


1.000001 


325 


1251 


140 


1111 


6.03 


1.17 


1.000001 


375 


1399 


140 


1259 


6.16 


1.14 


1.000000 


400 


1475 


142 


1333 


6.76 


1.26 


1.000000 



Table 4: Hartree Fock energy, centre of mass kinetic energy, mass and magnetic moment of 
A ++ where A' is 350 MeV, A is 100 MeV; b' is 0.60, b is 0.84 using nine by nine matrix. 



E HF (MeV) 


T CM (MeV) 


M(MeV) 


a 


r a „(fm) 


N 


1354 


130 


1224 


6.15 


1.22 


0.998476 



Figures (2) and (3) show the wave functions and the single particle potential for A' = 
350 MeV but A = 350 MeV and 100 MeV . As expected, there is no change in large r but an 
insignificant variation in small r. This is reflected in Ehf an d n ( table 2). 

On the other hand, A' variation ( figures 4 and 5) is quite significant and larger in large r. 
This is reflected in E HF and \x ( table 3). 

So we can conclude that for finite system calculations, the value of the confinement param- 
eter A' plays a more important role in the Richardson Potential. 
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1 2 3 4 5 

r in fm 



Figure 2: Variation of the wave functions with r for two different values of the asymptotic freedom parameter 
A ( using seven by seven matrices ); the upper ones are G(r) and the lower ones are F(r). Here the dotted curve 
is for A = A' = 350 MeVand the solid one is for A = 100, A' = 350 MeV. 




1 2 3 4 5 6 



r in fm 

Figure 3: Variation of the single particle potential with r for two different values of the asymptotic freedom 
parameter A ( using seven by seven matrices ). Here the dotted curve is for A = A' = 350 MeVand the solid 
one is for A = 100, A' = 350 MeV. 
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1 2 3 4 5 

r in fm 



Figure 4: Variation of the wave functions with r for two different values of the confinement parameter A' ( 
using seven by seven matrices ); the upper ones are G(r) and the lower ones are F(r) . Here the dotted curve 
is for A' = 375 MeV, A = 350 MeV and the solid one is for A' = 325 MeV, A = 350 MeV . 



8 - 



7 - 




1 2 3 4 5 6 

r in fm 

Figure 5: Variation of the single particle potential with r for two different values of the confinement parameter 
A' ( using seven by seven matrices ) . Here the dotted curve is for A' = 375 MeV, A = 350 MeV and the 
solid one is for A' = 325 MeV, A = 350 MeV . 
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4 Result with Q 

Q~ is the triple s-quark state. There is no qualitative change in the procedure except for for 
the quark mass; m s is taken here as 150 MeV. Results for A' = A = 350 MeV is shown in 
table [5]. 

Table 5: Variation of Hartree Fock energy, centre of mass kinetic energy, mass and magnetic 
moment of fl~ with oscillator parameter b where A = A' = 350 MeV, b' is 0.60. 



b 


E HF (MeV) 


7cm (MeV) 


M(MeV) 


tL 


r a i;(fm) 


N 


0.60 


1639 


176 


1463 


-1.95 


0.836879 


1.000001 


0.62 


1638 


175 


1463 


-1.95 


0.836778 


1.000000 


0.64 


1638 


174 


1464 


-1.95 


0.836879 


1.000000 


0.66 


1638 


172 


1466 


-1.95 


0.836854 


1.000001 


0.68 


1638 


171 


1467 


-1.95 


0.836882 


1.000000 


0.70 


1637 


169 


1468 


-1.95 


0.83643 


0.999999 


0.72 


1639 


167 


1472 


-1.95 


0.83695 


1.00001 



We improve the fit by taking different values for A' and A; specifically for A' = 350 MeV 
and A = 100 MeV. 

Table 6: Hartree Fock energy, centre of mass kinetic energy, mass and magnetic moment of VT 
where A' is 350 MeV, A is 100 MeV; b' is 0.60, b is 0.70. 



E HF (MeV) 


3cm (MeV) 


M(MeV) 


Ml 
MO 


r a „(fm) 


N 


1721 


165 


1556 


-1.92 


0.849015 


0.999999 



The centre of mass kinetic energy T cm is rather high resulting a lower mass of VL 
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5 Conclusions and summary 



In summary we have shown that a RHF calculation can be done for the energies and magnetic 
moments of the simplest baryons Q~ and A ++ . We employed the oscillator basis with good 
convergence. Optimized oscillator parameters are used for the large and small components of 
the RHF wave functions. We belive that this model may be good base for further investigations 
of baryon properties. The center of mass correction is around 150 MeV for these baryons (about 
10% of E HF ). Quark masses chosen are 4 and 150 MeV for u and s quark respectively. For q q 
interaction we chose Richardson potential. We separated out the confinement and asymptotic 
freedom scale parameter (A' and A ). From the best fit of the energies of Q~ and A ++ we find 
A' = 350 MeV and A = 100 MeV. We believe these are realistic values. Above parameters 
are now used in strange quark matter calculation ([21]). 
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